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We discuss the intriguing possibility that dark energy may change its equation of state in situations 
where large dark energy fluctuations are present. We show indications of this dynamical mutation 
in some generic models of dark energy. 



00 

o 

O 
(N 

a 



Oh' 

6 



> 
00 

m 
d 

o 



PACS numbers: 98.80.Cq 



INTRODUCTION 



There is aniple evidence that supports the existence of 
dark matter (DM) and dark energy (DE) in the universe 
Due to its charge neutrahty and dust-hke equation 
of state (i.e., neghgible pressure), dark matter starts to 
cluster gravitationally very early in the history of the 
Universe, and is crucial for the formation of large scale 
structure. Dark energy, on the other hand, becomes rele- 
vant only more recently, and is presumed to be a smooth 
component with a negative equation of state in order to 
fuel the accelerated expansion of the Universe . 

At the background level, dark energy (or any other 
fluid) is completely determined by its equation of state 
w = Pel Pe, where Pe is the pressure and pe is the energy 
density of dark energy. Already at this level dark energy 
can affect large scale structures [l] - see also [3| for the 
effect in different parameterizations of the equation of 
state of dark energy. 

However, if dark energy is in fact a manifestation of a 
dynamical mechanism such as a scalar field, then it will 
also develop inhomogeneities due to its gravitational in- 
teractions with itself and with dark matter [^. In linear 
perturbation theory, besides the energy density pertur- 
bation 5p, we need two extra degrees of freedom to char- 
acterize cosmological perturbations: the pressure pertur- 
bation 5p and the scalar anisotropic stress tt 0, Q ■ Alter- 
natively, one can also use the velocity potential 9 = V -v 
and the anisotropic stress [1]. 

The inhomogeneities of dark energy are often quite 
small, particularly in the case of ordinary (canonical) 
scalar field models with almost A-like behaviour - that is, 
when w c± — 1 5; In fact, as u" — > — 1 the perturba- 

tions in all dark energy models are suppressed in relation 
to those of dark matter. However, if that is not the case 
then the dark energy density contrast 5e = 5pe/ Pe can 
be either small or large, depending on the pressure per- 
turbations of dark energy [111 ] . 

Here we present further evidence of an intriguing 
possibility which was first pointed out in supergravity- 
motivated scalar field models of dark energy IJ, l3| : that 



dark energy can mutate into a fluid with clustering prop- 
erties similar to those of dark matter. We will show that 
this effect is a generic feature of dark energy, and that 
it has a simple origin: when pressure perturbations are 
large, the effective equation of state inside a collapsed 
region can be completely different from the equation of 
state of its homogeneous component. 



GRAVITATIONAL COLLAPSE WITH DARK 
ENERGY 



In the following, in order to specify the properties 
of the dark energy perturbations we will neglect the 
anisotropic stress. Moreover, we will characterize the 
pressure perturbation in a simplified manner, using the 
so-called effective (or non-adiabatic) sound velocity [14| . 
defined as c^^j = 5pe/5pe^ which we will assume to be a 
function of time only. Notice that this assumption lacks 
a formal basis in perturbation theory, since 8pe is a per- 
turbed variable whose time and spatial dependences can 
be, and often are, independent of the variations of 5pe- 
Nevertheless, in a particular gauge (the so-called "rest 



frame" of the fluid, where Tg = 0), the effective sound 

speed coincides with the phase velocity of linear relativis- 
ts 



tic perturbations, c\ Il5| 



Describing the pressure perturbation as 5pe — c^yj 5pe 
allows us to treat a wide variety of dark energy models 
and, crucially, it also allows us to compute non-linear 
structure formation using the Spherical Collapse model 
(SC) Furthermore, in this case the SC equations 

(derived in a simplified relativistic framework) are identi- 
cal to the equations of pseudo-Newtonian cosmology 17 1 . 
This means that the physics of gravitational collapse of 
structures such as galaxy clusters is well described within 
this framework. 

Consider then, in the spirit of the SC model, a spheri- 
cally symmetric region of constant dark energy overden- 
sity (the so-called "top-hat" density profile.) Let us call 
Pg = Pe+5pe andpg = Pe+5pe the energy density and the 
pressure of this region, which are modified with respect 
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to the corresponding background quantitites, pe and Pe 
by the perturbations 5pe and 5pe- The equation of state 
w is defined as the ratio of the total pressure to the to- 
tal energy density, and hence it will be different for the 
background and the interior of the collapsed region. A 
simple calculation shows that the equation of state inside 
the collapsed region, w'^, is given by: 
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For small density contrasts \Se\ ^ 1, the equation of 
state inside the overdense region does not change ap- 
preciably. However, if c^^y ^ w, then in the nonlinear 
regime, where S > 1 (halos), there could be a substan- 
tial modification in w'^ with respect to the background 
equation of state. Even in underdense regions (voids), 
where 5^—1, there could be large modifications of the 
equation of state. Hence, in principle dark energy could 
even effectively mutate into dark matter inside halos and 
voids. 

The above argument is completely general. What re- 
mains to be shown is whether there are models in which 
this dramatic situation is actually realized. This requires 
a non-linear analysis of the evolution of pertubations for 
two gravitationally coupled fluids, dark energy and dark 
matter (we will neglect radiation and baryons in what 
follows). Unfortunately, at present there are no totally 
rigorous methods for performing this analysis - except 
in the case of canonical scalar fields, but even then only 
approximately (ol. [loj. 

Here we employ a generalization of the SC model for 
the case of a relativistic fluid with pressure. In the next 
section we present the relevant equations and mention 
under which conditions they are equivalent to a pseudo- 
Newtonian approach. We will then analyse the evolution 
of the coupled system of perturbations for a wide variety 
of dark energy models for which the pressure perturba- 
tions are characterized by some homogeneous effective 
sound speed. 



NON-LINEAR EVOLUTION 

We define H = d/a and h = rjr as the expansion 
rates for the background (a is the scale factor) and for 
the perturbed region (r is the size of the collapsing re- 
gion), respectively. In what follows we work within two 
assumptions, namely, that there is no non-gravitational 
interaction between DE and DM, and that the total en- 
ergy of both DE and DM contained in the collapsed re- 
gion is constant. The possibility of including DM-DE in- 
teractions in the study of structure formation was studied 
m and a discussion of possible outflow of energy 



from the collapsed region can be found in [IJ, [2 



for the perturbed region for a fluid species j: 

^ = -3i/(l + w;,) , 4 = -3Ml + ^«p, (2) 

we obtain, for the density contrasts: 

6, = -3(1 + 5,) [h{\ + - H(l + ,«,)] . (3) 

Obviously, for matter we have Wm = w'^ — 0. 

Using Eq. ([T]) and that the local expansion rate is 
related to the velocity field in the perturbed region by 
h ^ H + 9 /3a, we arrive at: 

Q 

5j+m{c%jf-wj)5j + - [(1 + w,) + + c^j,}f)5j\ = . 

(4) 

The equation that determines the evolution of 9 comes 
from the "acceleration" in the perturbed region: 



-H9 + — + -H^a[Qm5n 
6a I 



+ (1 + 30^^^)^!^] =0. (5) 



Notice that there is only one equation for the peculiar 
velocity, even in the case of 2 fluids. This is clearly nec- 
essary, because in the SC model it is a single spherically 
symmetric region that detaches from the background, 
with a peculiar expansion rate given by /i = -I- 9 /ia. 
When there is pressure and pressure gradients, relativis- 
tic corrections almost surely break this identity between 
the velocities of the fluids, which means that the SC 
model with a top-hat profile is inconsistent with a fully 
relativistic calculation. The assumption that c^j^ is only 
time dependent implies that 5p is also only time depen- 
dent - which ultimately guarantees the validity of the top 
hat SC model. 

The SC equations capture many features of the gravi- 
tational physics one expects to find at the scale of the col- 
lapsed structures we see today. This is because the exact 
same equations can be derived from a pseudo-Newtonian 
treatment of perturbations when gradients of pressure 
can be neglected - as is the present case, of a top hat 
profile. When linearized, the ensuing equations corre- 
spond to the linear equations from General Relativity in 



the case 



— for sub-horizon scales 11? 



Further- 



more, even when c^y^ ^ 0, although the equations are 
not equivalent anymore, still the growing modes of the 
linearized pseudo-Newtonian perturbations are identical 
to the growing modes of the linearized relativistic per- 
turbations [17|. 



Using the continuity equations for the background and 



NUMERICAL SOLUTIONS 

We solved the coupled differential equations for Jg, 5m 
and 6*, for a few representative models of dark energy. In 
our approximation, a dark energy model is determined by 
its background equation of state and the effective speed of 
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FIG. 1: Clustered dark energy equation of state for different 
cases of effective speed of sound for background w — —0.8: 
clff = w (solid line), C^ff = 1 (dashed line), c^ff — (dot- 
dashed line), ciff~ —w (dotted line) and c^jj = —1 (double 
dot-dashed line). The instant of turnaround (/i = 0) is z = 0.3 
for c^j j- = w, z — 0.5 for both c^jf = —w and c^ff — 1, and 
z — 0.6 for c^ff = (for c^^j = —1 there is no turnaround.) 



sound. In a previous paper we investigated the particular 
case where the sound speed of dark energy is equal to its 
equation of state, c^j^ — w, and therefore w'^ = w. In 
that case, one can see clearly from Eq. ([T]) that there is 
no mutation 



ill 



In this letter we expand our previous analysis to the 



following cases: c; 



I.e. 



eff 



0,c; 



eff 



-land clff 



—w. The first case is motivated by the common situation 
when dark energy is modelled by a canonical scalar field, 
since in the gauge corresponding to the rest frame of the 
scalar field the effective sound speed is Cg^j = 
The second case (null pressure perturbations) can occur 
in so-called "silent quartessence" models [2l'|. The third 
case represents a perturbation with behaviour close to a 
cosmological constant and the last case reproduces the so- 
called generalized Chaplygin gas in a certain limit [a = 

1) di. 

We use adiabatic initial conditions for the perturba- 
tions [SI = (1 -I- w)6l^] at a redshift z — 1000, an initial 
velocity field coincident with the Hubble flow (6'' = 0), 
Ha = 72 km s^^ Mpc"\ Qm = 0.25 and = 0.75. 
In the examples shown below we fixed the background 
equation of state of dark energy at w = —0.8. 

In Fig. 1 we show the values of the clustered equation 
of state Wc as a function of the redshift. The initial con- 
dition Sl^ was chosen such that the dark energy perturba- 
tion Se ~ today. The perturbation in dark matter 
is typically one to two orders of magnitude larger, which 
is consistent with the typical density contrast in galaxy 
clusters, for which Sm ^ O(IO^). 

As discussed above, for the case c^j^ = w we obtain no 
mutation in the equation of state in the perturbed region. 



FIG. 2: Clustered dark energy equation of state for different 
cases of effective speed of sound for background w = —0.99. 
Lines are the same as in Figure 1. The cases of c^^j — w and 



1, as well as c^ft — 1 and c^tt = —w, lie on top of 



each other. Turnaround occurs at approximately z ~ 0.6 in 
all cases. 



However, substantial modifications are found for other 
possibilities of c^ff- The largest modifications arise for 
Cgyj — 1, which in our example is very similar to c^^j = 
—w. In these two cases, even a complete metamorphosis 
of dark energy into a fluid which clusters as strongly as 
dark matter is possible at recent epochs, due to the large 
perturbations inside the collapsed region. The behaviour 
for Cgjj = — 1 is also easily understood from Eq. ([T]), 
since this is the only case where Cg^j — w < 1. 

As expected, the effect of mutation is greatly reduced 
for a background equation of state close to that of a cos- 
mological constant, independent of c^ff- We illustrate 
that fact in Fig. 2, where we show the equation of state 
inside the collapsing region for different models of clus- 
tered dark energy, in the case of a background equation 
of state w = —0.99, and with the same initial conditions 
that were used in Fig. 1. Nevertheless, even in this case 
large density contrasts can still arise in the dark energy 
component, which lead to the mutation of dark energy. 



CONCLUSIONS 

We have shown that it is possible to change radically 
the clustering properties of dark energy in collapsed re- 
gions (halos and voids.) We exemplifled this behaviour 
with a few models for the dark energy perturbations, and 
showed that it happens not only in scalar field models, 
but also in generic models of dark energy - in particular 
the Generalized Chaplygin Gas and Silent Quartessence 
models. 

Since the physics of most observed collapsed struc- 
tures, such as galaxy clusters, is well approximated 



4 



by quasi-Newtonian physics, this dynamical mutation 
should be a general phenomenom. Clearly, this is a cru- 
cial issue for all attempts to compute the influence of 
dark energy on the formation of large scale structures. 
More detailed studies, including a relativistic approach 
and using different realistic parameterizations of the dark 
energy equation of state are currently under way 17| . 



Acknowledgments 

We would like to thank loav Waga for many fruitful 
discussions. This work has been supported by FAPESP 
grants 04/13668-0 (L.R.A. and R.R.) and 05/00554-0 
(R.C.B.), a CNPq grant 309158/2006-0 (R.R.) and a 
CAPES grant (L.L.). 



[1] See, e.g., U. Seljak, A. Slozar and P. McDonald, JCAP 

0610, 014 (2006). 
[2] For reviews, see e.g. T. Padmanabhan, Phys. Kept. 380: 

325 (2003); P. J. E. Peebles and B. Ratra, Rev. Mod. 

Phys. 75: 559 (2003); V. Sahni and A. A. Statobinsky, 

Int. J. Mod. Phys. D 15, 2105 (2006). 
[3] See, e.g., L. M. Wang and P. J. Steinhardt, Ap. J. 508: 

483 (1998); R. A. Battye and J. Weller, Phys. Rev. D68; 

083506 (2003); E. V. Linder, Phys. Rev. D72: 043529 

(2005). 



[4] L. Liberato and R. Rosenfeld, JCAP 0607, 009 (2006). 
[5] K. Coble, S. Dodelson and J. A. Frieman, Phys. Rev. 

D55: 1851 (1997). 
[6] J. M. Bardeen, Phys. Rev. D 22, 1882 (1980). 
[7] H. Kodama and M. Sasaki, Progr. Theor. Phys. Suppl. 

78: 1 (1984). 

[8] C.-P. Ma and E. Bertschinger, Ap. J. 455: 7 (1995). 
[9] S. Dutta and I. Maor, Phys. Rev. D 75, 063507 (2007). 
[10] D. F. Mota, D. J. Shaw and J. Silk, [arXiv: 0709 . 2227,1 . 
[11] L. R. Abramo, R. C. Batista, L. Liberato and R. Rosen- 
feld, JCAP 0711:12 (2007). 
[12] D. F. Mota and C. van de Bruck, Astron. Astrophys. 
421:71 (2004). 

[13] N. J. Nunes and D. F. Mota, Mon. Not. R. Ast. Soc. 368, 
751 (2006) 

[14] W. Hu, Astrophys. J. 506, 485 (1998). 
[15] J. Garriga and V. Mukhanov, Phys. Lett. B 458, 219 
(1999). 

[16] J. E. Gunn and J. R. Gott III, Ap. J. 176: 1 (1972). 

[17] L. R. Abramo, R. C. Batista, L. Liberato and R. Rosen- 
feld, in preparation. 

[18] R. R. R. Reis, Phys. Rev. D 67, 087301 (2003), Erratum- 
ibid. D 68, 089901 (2003). 

[19] M. Manera and D. F. Mota, Mon. Not. R. Ast. Soc. 371, 
1373 (2006). 

[20] N. J. Nunes, A. C. da Silva and N. Aghanim, Astron. 

Astrophys. 450:899 (2006). 
[21] L. Amendola, I. Waga and F. FineUi, JCAP 0511:009 

(2005). 

[22] L. Amendola, F. Fineli, C. Burigana and D. Carturan, 
JCAP 0307:005 (2003). 



